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Predicting the cosmological constant with the scale-factor cutoff measure 
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It is well known that anthropic selection from a landscape with a flat prior distribution of cosmo- 
logical constant A gives a reasonable fit to observation. However, a realistic model of the multiverse 
has a physical volume that diverges with time, and the predicted distribution of A depends on how 
the spacetime volume is regulated. We study a simple model of the multiverse with probabilities reg- 
ulated by a scale-factor cutoff, and calculate the resulting distribution, considering both positive and 
negative values of A. The results are in good agreement with observation. In particular, the scale- 
factor cutoff strongly suppresses the probability for values of A that are more than about ten times 
the observed value. We also discuss several qualitative features of the scale-factor cutoff, including 
aspects of the distributions of the curvature parameter fi and the primordial density contrast Q. 



I. INTRODUCTION 

The present understanding of inflationary cosmology 
suggests that our universe is one among an infinite num- 
ber of "pockets" in an eternally inflating multiverse. 
Each of these pockets contains an infinite, nearly homoge- 
neous and isotropic universe and, when the fundamental 
theory admits a landscape of metastable vacua, each may 
be characterized by different physical parameters, or even 
different particles and interactions, than those observed 
within our pocket. Predicting what physics we should 
expect to observe within our region of such a multiverse 
is a major challenge for theoreticalphysics. (For recent 
reviews of this issue, see e.g. [3, 0,0a HQ-) 

The attempt to build a calculus for such predictions is 
complicated in part by the need to regulate the diverg- 
ing spacetime volume of the multiverse. A number of 
different approaches to this measure problem has been 
explored: a cutoff at a fixed global time [aullgl 1 , the so- 
called "gauge-invariant" measures [Til . [l2| , where differ- 
ent cutoff times are used in different pockets in order to 
make the measure approximately time-par ametrization 
invariant, the pocket-based measure [HI, E3, EH, EH, 
which avoids reference to global time by focusing on 
pocket abundances and regulates the diverging volume 
within each pocket with a spherical volume cutoff, and 
finally the causal patch measures [13, EH , which restrict 
consideration to the spacetime volume accessible to a sin- 
gle observer. 2 Different measures make different obser- 
vational predictions. In order to decide which, if any, is 



Much of the early work sought to calculate the relative vol- 
umes occupied by different pockets on hypersurfaces of constant 
time [9, 10]. In Ref. [f| the probabilities were expressed in terms 
of the fluxes appearing in the Fokker-Planck equation for eternal 
inflation. In most (but not all) cases, this method is equivalent 
to imposing a cutoff at a constant time. The prescription of a 
global time cutoff was first explicitly formulated in @| . 
We also note the recent measure proposals in Refs. |l9t |2(J . 
Observational predictions of these measures have not yet been 
worked out, so we shall not discuss them any further. 



on the right track, one can take an empirical approach, 
working out the predictions of candidate measures and 
comparing them with the data. In this spirit, we investi- 
gate one of the simplest global-time measure proposals: 
the scale-factor cutoff measure. 

The main focus of this paper is on the prediction of 
the cosmological constant A [2l], [H, [H, 0, H3, , which 
is arguably a major success of the multiverse picture. 
Most calculations of the distribution of A in the litera- 
ture [H [H [H, [13, [H, Hi] do not explicitly specify the 
measure, but in fact correspond to using the pocket-based 
measure. The distribution of positive A in a causal-patch 
measure has also been considered (2fj- The authors of 
Ref. [3fj| emphasize that the causal-patch measure gives 
a strong suppression for values of A more than about 
ten times the observed value, while anthropic constraints 
alone might easily allow values 1000 times larger than 
observed, depending on assumptions. Here, we calculate 
the distribution for A in the scale-factor cutoff measure, 
considering both positive and negative values of A, and 
compare our results with those of other approaches. We 
find that our distribution is in a good agreement with 
the observed value of A, and that the scale-factor cutoff 
gives a suppression for large positive values of A that is 
very similar to that of the causal-patch measure. 

We also show that the scale-factor cutoff measure is 
not afflicted with some of the serious problems arising 
in other approaches. For example, another member of 
the global time measure family — the proper-time cut- 
off measure — predicts a population of observers that 
is extremely youth-dominated [3J 0, [Hj] . Observers who 
take a little less time to evolve are hugely more numerous 
than their slower-evolving counterparts, suggesting that 
we should most likely have evolved at a very early cos- 
mic time, when the conditions for life were rather hostile. 
This counter-factual prediction is known as "the young- 
ness paradox" . Furthermore, the "gauge-invariant" and 
pocket-based measures suffer from a U Q catastrophe," 
exponentially preferring either very large or very small 
values of the primordial density contrast Q [33L l34j . In 
fact, this problem is not restricted to Q — there are simi- 
lar expectations for the gravitational constant G [3|| . We 
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show that the youngness bias is very mild in the scale- 
factor cutoff, and that there is no Q (or G) catastrophe. 
We also describe qualitative expectations for the distri- 
butions of Q and of the curvature parameter Q. 

This paper is organized as follows. In section |TT] we 
describe the scale-factor cutoff, commenting on its more 
salient features including its very mild youngness bias 
and aspects of the distributions of Q and O. In sec- 
tion IIIII we compute the probability distribution of A, 
calculating it first for the pocket-based measure, repro- 
ducing previous results, and then calculating it for the 
scale-factor cutoff. In both cases we study positive and 
negative values of A. Our main results are summarized 
in section IIVI Finally, we include two appendices. In 
appendix |A] we consider the possibility that the land- 
scape splits into several disconnected sectors, and show 
that even in this situation the scale-factor cutoff measure 
is essentially independent of the initial state of the uni- 
verse. Appendix [B] contains an analysis of the evolution 
of the collapse density threshold, along with a description 
of the linear growth function of density perturbations. 



II. THE SCALE-FACTOR CUTOFF 

A. Global time cutoffs 

To introduce a global time cutoff, we start with a patch 
of a spacelike hypersurface E somewhere in the inflat- 
ing part of spacetime, and follow its evolution along the 
congruence of geodesies orthogonal to E. The spacetime 
region covered by this congruence will typically have infi- 
nite spacetime volume, and will include an infinite num- 
ber of pockets. In the global-time cutoff approach we 
introduce a time coordinate t, and restrict our attention 
to the finite spacetime region r(E,i c ) swept out by the 
geodesies prior to t = t c , where t c is a cutoff which is 
taken to infinity at the end of the calculation. The rel- 
ative probability of any two types of events A and B is 
then defined to be 



p(A) _ n(A,r(£,f c )) 

t c — ^OC 



p(B) 



i(B,r(E,t c )) ' 



(1) 



where n(A,T) and n(_B,r) are the number of events of 
types A and B respectively in the spacetime region T. 
In particular, the probability Pj of measuring parame- 
ter values corresponding to a pocket of type j is pro- 
portional to the number of independent measurements 
made in that type of pocket, within the spacetime region 
r(E, t c ), in the limit t c — > oo. 

The time coordinate t is "global" in the sense that 
constant-time surfaces cross many different pockets. 
Note however that it does not have to be global for the 
entire spacetime, so the initial surface E does not have 
to be a Cauchy surface for the multiverse. It need not 
be monotonic, either, where for nonmonotonic t we limit 
r(E,i c ) to points along the geodesies prior to the first 
occurrence of t = t r . 



As we will discuss in more detail in appendix [Al prob- 
ability distributions obtained from this kind of measure 
are independent of the choice of the hypersurface E. 3 
They do depend, however, on how one defines the time 
parameter t. To understand this sensitivity to the choice 
of cutoff, note that the eternally inflating universe is 
rapidly expanding, such that at any time most of the 
volume is in pockets that have just formed. These pock- 
ets are therefore very near the cutoff surface at t = t c , 
which explains why distributions depend on exactly how 
that surface is drawn. 

A natural choice of the time coordinate t is the proper 
time t along the geodesic congruence. But as we have 
already mentioned, and will discuss in more detail in 
the following subsection, this choice is plagued with the 
youngness paradox, and therefore does not yield a sat- 
isfactory measure. Another natural option is to use the 
expansion factor a along the geodesies as a measure of 
time. The scale-factor time is then defined as 



t = In a . 



(2) 



The use of this time parameter for calculating probabil- 
ities is advocated in Ref. [36j and is studied in various 
contexts in Refs. [lj|, [gj, and [j|. 4 It amounts to 
measuring time in units of the local Hubble time JT _1 , 



dt = Hdr . 



(3) 



The scale-factor cutoff is imposed at a fixed value of t = 
i c , or, equivalently, at a fixed expansion factor a c . 

The term "scale factor" is often used in the context of 
homogeneous and isotropic spaces, but it is easily gener- 
alized to spacetimes with no such symmetry. In the gen- 
eral case, the scale-factor time can be defined by Eq. ([3]) 
with 



ff=(l/3)<. 



(4) 



where is the four- velocity vector along the 

geodesies. This definition has a simple geometric mean- 
ing, which can be seen by imagining that the congruence 
of geodesies describes the flow of a "dust" of test parti- 
cles. If the dust of particles is assumed to have a uniform 



3 Here, and in most of the paper, we assume an irreducible land- 
scape, where any metastable inflating vacuum is accessible from 
any other such vacuum through a sequence of transitions. Alter- 
natively, if the landscape splits into several disconnected sectors, 
each sector will be characterized by an independent probability 
distribution and our discussion will still be applicable to any of 
these sectors. The distribution in case of a reducible landscape 
is discussed in appendix lAl 

4 The measure studied in Ref. [36l | is a comoving-volume measure 
on surfaces of constant scale-factor time; it is different from the 
scale-factor cutoff measure being discussed here. In particular, 
the former measure has a strong dependence on the initial state 
at the hypersurface S. Our measure is very similar to one stud- 
ied in Ref. @], which is called the "pseudo-comoving volume- 
weighted measure." 
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density po on the initial surface E, then the four-current 
of the dust can be described by j^{x) — p{x)u^{x), where 
p = po on S. Conservation of the current then implies 
that u^d^p + = 0, which with Eqs. © and 

implies that 

D T In p = = — 3D T t , (5) 

where D T = u^d^ is the derivative with respect to 
proper time along the geodesies. The solution is then 
p = poe~ 3t . From Eq. ((2|) we then have a cx p -1 ^ 3 , so 
the scale-factor cutoff is triggered when the density p(x) 
of the dust in its own rest frame drops below a certain 
specified level. 

The divergence of geodesies during inflation or homo- 
geneous expansion can be followed by convergence during 
structure formation or in regions dominated by a negative 
cosmological constant. The scale- factor time then ceases 
to be a good time variable, but this does not preclude 
one from using it to impose a cutoff. A geodesic is termi- 
nated when the scale factor first reaches the cutoff value 
a c . If the scale factor turns around and starts decreasing 
before reaching that value, we continue the geodesic all 
the way to the crunch. When geodesies cross we can still 
define the scale factor time along each geodesic according 
to Eqs. ^ and then one includes a point in r(E, t c ) 
if it lies on any geodesic prior to the first occurrence of 
t = t c on that geodesic. 

To facilitate further discussion, it will be useful to re- 
view some general features of eternally inflating space- 
times, and how they are reflected in proper time and 
scale-factor time slicings. Regions of an eternally inflat- 
ing multiverse may evolve in two distinct ways. In the 
case of quantum diffusion [3?], HH , inflation is driven by 
the potential energy of some light scalar fields, the evo- 
lution of which is dominated by quantum fluctuations 
and is described by the Fokker-Planck equation (see e.g. 
Ref. [HI). Pockets form when the scalar field(s) fluctuate 
into a region of parameter space where classical evolution 
dominates, and slow-roll inflation ensues. One can de- 
fine spacelike hypersurfaces separating t he q uantum and 
classical regimes (see for example Ref. [HI]); which we 
denote by E g . In universes like ours, slow-roll inflation is 
followed by thermalization (reheating) and the standard 
post-inflationary evolution. We denote the hypersurface 
of thermalization, which separates the inflationary and 
post-inflationary epochs, as E*. 

The multiverse may also (or instead) feature massive 
fields associated with large false- vacuum energies. Evolu- 
tion is then governed by bubble nucleation through quan- 
tum tunneling [39l . |40T ] and can be described to good ap- 
proximation by a suitable master equation [ill . fl5| . The 
tunneling may proceed into another local minimum, into 
a region of quantum diffusion, or into a region of classical 
slow-roll inflation. In the latter case, the bubble interiors 
have the geometry of open FRW universes [13]. Bubbles 
of interest to us here have a period of slow-roll inflation 
followed by thermalization. The role of the hypersurface 
E g is played in this case by the surface separating the ini- 



tial curvature-dominated regime and the slow-roll regime 
inside the bubble. The differences between quantum dif- 
fusion and tunneling are not important for most of the 
discussion below, so we shall use notation and terminol- 
ogy interchangeably. 

The number of objects of any type that have formed 
prior to some time t is proportional to e 7 *, where 7 is the 
largest eigenvalue of the physical-volume Fokker-Planck 
or master equation. This is because the asymptotic be- 
havior is determined by the eigenstate with the largest 
eigenvalue. Similarly, the physical volume that thermal- 
izes into pockets of type j between times t and t + dt has 
the form 

dV^ =C J e lt dt, (6) 

where Cj is a constant that depends on the type of 
pocket. (This was derived in Ref. @ for models with 
quantum diffusion and in Refs. [43j] and [Hj] for models 
with bubble nucleation.) 

The value of 7 in Eq. ^ is the same for all pockets, 
but it depends on the choice of time variable t. With a 
proper-time slicing, it is given by 

7 ~ 3ff max (* = t) , (7) 

where H max is the expansion rate of the highest-energy 
vacuum in the landscape, and corrections associated with 
decay rates and upward tunneling rates have been ig- 
nored. In this case the overall expansion of the multiverse 
is driven by this fastest-expanding vacuum, which then 
"trickles down" to all of the other vacua. With scale- 
factor slicing, all regions would expand as a 3 = e 3t if it 
were not for the continuous loss of volume to terminal 
vacua with negative or zero A. Because of this loss, the 
value of 7 is slightly smaller than 3, and the difference 
is determined mostly by the rate of decay of the slowest- 
decaying (dominant) vacuum in the landscape (44|, 

7«3-kz3 (f = lna). (8) 

Here. 

n D = (4tt/3) Y d /H% , (9) 

where Yd is the decay rate of the dominant vacuum per 
unit spacetime volume, and Hd is its expansion rate. The 
vacuum decay rate is typically exponentially suppressed, 
so for the slowest-decaying vacuum we expect it to be 
extremely small. Hence, 

3-7 <1. (10) 

B. The youngness bias 

As we have already mentioned, the proper-time cutoff 
measure leads to rather bizarre pre dictions, collectively 
known as the youngness paradox |3ll . 0, l32| . With proper 
time slicing, Eqs. ^ and ([7]) tell us that the growth of 
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volume in regions of all types is extremely fast, so at 
any time the thermalized volume is exponentially domi- 
nated by regions that have just thermalized. With this 
super-fast expansion, observers who take a little less time 
to evolve are rewarded by a huge volume factor. This 
means most observers form closer to the cutoff, when 
there is much more volume available. Assuming that 
H max is comparable to Planck scale, as one might ex- 
pect in the string theory landscape, then observers who 
evolved faster than us by At = 10 9 years would have 
an available thermalized volume which is larger than the 
volume available to us by a factor of 



,3H n 



exp(10 6 



(11) 



Unless the probability of life evolving so fast is sup- 
pressed by a factor greater than exp(10 60 ), then these 
rapidly evolving observers would outnumber us by a huge 
factor. Since these observers would measure the cos- 
mic microwave background (CMB) temperature to be 
T = 2.9 K, it would be hard to explain why we measure 
it to be T = 2.73 K. Note that because H max AT appears 
in the exponent, the situation is qualitatively unchanged 
by considering much smaller values of i? max or At. 

The situation with a scale-factor cutoff is very differ- 
ent. To illustrate methods used throughout this paper, 
let us be more precise. Let At denote the interval in 
scale- factor time between the time of thermalization, t*, 
and the time when some class of observers measures the 
CMB temperature. A time cutoff excludes the count- 
ing of observers who measure the CMB temperature at 
times later than t c , so the number of counted observers 
is proportional to the volume that thermalizes at time 
t* < t c — At. (For simplicity we focus on pockets that 
have the same low-energy physics as ours.) The volume 
of regions thermalized per unit time is given by Eq. ([6]). 
During the time interval At, some of this volume may de- 
cay by tunneling transitions to other vacua. This effect 
is negligible, and we henceforth ignore it. For a given At, 
the thermalized volume available for observers to evolve, 
as counted by the scale-factor cutoff measure, is 



V(At) 



i c -At 



(12) 



To compare with the results above, consider the rela- 
tive amounts of volume available for the evolution of two 
different civilizations, which form at two different time 
intervals since thermalization, At\ and At 2 : 



V(Ah 
V{At 2 



3 T(A* a -Atx) = (a2 / ai) 7 ; 



(13) 



where aj is the scale factor at time + Atj. Thus, taking 
7 ~ 3, the relative volumes available for observers who 
measure the CMB at the present value (T — 2.73 K), 
compared to observers who measure it at the value of 
10 9 years ago (T = 2.9 K), is given by 



V(2.73 K) /2.73K 



V(2.9 K) 



2.9 K 



0.8. 



(14) 



Thus, the youngness bias is very mild in the scale-factor 
cutoff measure. Yet, as we shall see, it can have interest- 
ing observational implications. 



C. Expectations for the density contrast Q and the 
curvature parameter fl 

Pocket-based measures, as well as "gauge-invariant" 
measures, suffer from a U Q catastrophe" where one ex- 
pects to measure extreme values of the primordial density 
contrast Q. To see this, note that these measures expo- 
nentially prefer parameter values that generate a large 
number of e-folds of inflation. This by itself does not ap- 
pear to be a problem, but Q is related to parameters that 
determine the number of e-folds. The result of this is a 
selection effect that exponentially prefers the observation 
of either very large or very small values of Q, depending 
on the model of inflation and on which inflationary pa- 
rameters scan (i.e., which parameters vary significantly 
across the landscape) (3^, Sj. On the other hand, we ob- 
serve Q to lie comfortably in the middle of the anthropic 
range [45|], indicating that no such strong selection ef- 
fect is at work. 5 Note that a similar story applies to the 
magnitude of the gravitational constant G [351 ] . 

With the scale-factor cutoff, on the other hand, this is 
not a problem. To see this, consider a landscape in which 
the only parameter that scans is the number of e-folds of 
inflation; all low-energy physics is exactly as in our uni- 
verse. Consider first the portions of the hypersurfaces S g 
that begin slow-roll inflation at time t q in the interval dt q . 
These regions begin with a physical volume proportional 
to e 7 * 9 dt q , and those that do not decay grow by a factor 
of e 3N " before they thermalize at time t* — t q + N e . If 
kj is the transition rate out of the slow-roll inflationary 
phase (as defined in Eq. ([9])), then the fraction of volume 
that does not undergo decay is e~ K ' N,! . 

After thermalization at time t*, the evolution is the 
same in all thermalized regions. Therefore we ignore this 
common evolution and consider the number of observers 
measuring a given value of N e to be proportional to the 
volume of thermalization hypersurfaces that appear at 
times earlier than the cutoff at scale-factor time t c . This 
cutoff requires t* = t„ + N e < t c . Summing over all times 



tq gives 



P(N e ) cx e* 3 -"'^' 



-Ne 



e ltq dtq cx ( 



(3-7-K7)Ar e 



(15) 



Even though the dependence on N e is exponential, the 
factor 



7 - Ki PS K D - Kl 



(16) 



Possible resolutions to this problem have been proposed in 

Rets, mm be si- 
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is exponentially suppressed. Thus we find P(N e ) is a very 
weak function of N e , and there is not a strong selection 
effect for a large number of e-folds of slow-roll inflation. 
In fact, since the dominant vacuum D is by definition 
the slowest-decaying vacuum, we have hi > up. Thus 
the scale-factor cutoff introduces a very weak selection 
for smaller values of N e . 6 

Because of the very mild dependence on N e , we do 
not expect the scale-factor measure to impose significant 
cosmological selection on the scanning of any inflationary 
parameters. Thus, there is no Q catastrophe — nor is 
there the related problem for G — and the distribution 
of Q is essentially its distribution over the states in the 
landscape, modulated by inflationary dynamics and any 
anthropic selection effects. 

The distribution P(N e ) is also important for the ex- 
pected value of the curvature parameter Q. This is be- 
cause the deviation of Q from unity decreases during an 
inflationary era, 
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-2N. 



(17) 



Hence pocket-based and "gauge-invariant" measures, 
which exponentially favor large values of N e , predict a 
universe with extremely close to unity. The distribu- 
tions of f2 from a variety of models have been calculated 
using a pocket-based measure in Refs. [Hj] and [43| . 

On the other hand, as we have just described, the scale- 
factor cutoff measure does not significantly select for any 
value of N e . There will still be some prior distribution of 
N e , related to the distributions of inflationary parameters 
over the states in the landscape, but it is not necessary 
that N e be driven strongly toward large values (in fact, it 
has been argued that small values should be preferred in 
the string landscape, see e.g. Ref. [II]). Thus, it appears 
that the scale-factor cutoff allows for the possibility of 
a detectable negative curvature. The probability distri- 
bution of £1 in this type of measure has been discussed 
qualitatively in Ref. [48|; a more detailed quantitative 
analysis will be given elsewhere [49| . 



III. THE DISTRIBUTION OF A 

A. Model assumptions 

We now consider a landscape of vacua with the same 
low-energy physics as we observe, except for an essen- 
tially continuous distribution of possible values of A. Ac- 
cording to Eq. ([HI), the volume that thermalizes between 
times i* and t* + dt* with values of cosmological constant 



We are grateful to Ben Freivogel for pointing out to us the need 
to account for vacuum decay during slow-roll inflation. He has 
also suggested that this effect will lead to preference for smaller 
values of N e . 



between A and A + dA is given by 

dV* (A) = C{A)dAe lt "dU 



(18) 



The factor of C (A) plays the role of the "prior" distribu- 
tion of A; it depends on the spectrum of possible values 
of A in the landscape and on the dynamics of eternal 
inflation. The standard argument [U, [2H suggests that 
C(A) is well approximated by 



C(A) w const , 



(19) 



because anthropic selection restricts A to values that are 
very small compared to its expected range of variation in 
the landscape. The conditions of validity of this heuristic 
argument have been studied in simple landscape mod- 
els Q H3, [EH , with the conclusion that it does in fact 
apply to a wide class of models. Here, we shall assume 
that Eq. HU) is valid. 

Anthropic selection effects are usually characterized 
by the fraction of matter that has clustered in galax- 
ies. The idea here is that a certain average number of 
stars is formed per unit galactic mass and a certain num- 
ber of observers per star, and that these numbers are 
not strongly affected by the value of A. Furthermore, 
the standard approach is to assume that some minimum 
halo mass Mg is necessary to drive efficient star forma- 
tion and heavy element retention. Since we regulate the 
volume of the multiverse using a time cutoff, it is impor- 
tant for us to also track at what time observers arise. We 
assume that after halo collapse, some fixed proper time 
lapse At is required to allow for stellar, planetary, and 
biological evolution before an observer can measure A. 
Then the number of observers measuring A before some 
time r in a thermalizcd volume of size V* is roughly 



jV oc F(Mc, t — At)T4 



(20) 



where F is the collapse fraction, measuring the fraction 
of matter that clusters into objects of mass greater than 
or equal to Mq, at time t — At. 

Anthropic selection for structure formation ensures 
that within each relevant pocket matter dominates the 
energy density before A does. Thus, all thermalized re- 
gions evolve in the same way until well into the era of 
matter domination. To draw upon this common evolu- 
tion, within each pocket we define proper time r with 
respect to a fixed time of thermalization, t*. It is conve- 
nient to also define a reference time r m such that r m is 
much larger than the time of matter-radiation equality 
and much less than the time of matter-A equality. Then 
evolution before time r m is the same in every pocket, 
while after r rn the scale factor evolves as 



~a(r)- 
Here we have defined 




sinh 2/3 (§#at) forA>0 



,2/3/3 



Hat) 



for A < . 



H A = \p, 



(21) 



(22) 



6 



and use units with G = c = 1. The prefactors H~^ 2 ^ 
ensure that early evolution is identical in all thermalized 
regions. This means the global scale factor a is related 
to a by some factor that depends on the scale-factor time 

at which the region of interest thermalized. 

In the case A > 0, the rate at which halos accrete mat- 
ter decreases with time and halos may settle into galaxies 
that permit quiescent stellar systems such as ours. The 
situation with A < is quite different. At early times, 
the evolution of overdensities is the same; but when the 
proper time reaches r turn = tt/3H Ai the scale factor be- 
gins to decrease and halos begin to accrete matter at a 
rate that increases with time. Such rapid accretion may 
prevent galaxies from settling into stable configurations, 
which in turn would cause planetary systems to undergo 
more frequent close encounters with passing stars. This 
effect might become significant even before turnaround, 
since our present environment benefits from positive A 
slowing the collision rate of the Milky Way with other 
systems. 

For this reason, we use Eq. ([20]) to estimate the number 
of observers if A > 0, but for A < we consider two 
alternative anthropic hypotheses: 

A. we use Eq. (f^Uj) . but of course taking account of the 
fact that the proper time r cannot exceed r crunc h = 
2tt/3H a : or 

B. we use Eq. (f^Uj) . but with the hypothesis that the 
proper time r is capped at Tturn = tt/3Ha. 

Here r crunc h refers to the proper time at which a ther- 
malized region in a collapsing pocket reaches its future 
singularity, which we refer to as its "crunch." Anthropic 
hypothesis A corresponds to the assumption that life can 
form in any sufficiently massive collapsed halo, while an- 
thropic hypothesis B reflects the assumption that the 
probability for the formation of life becomes negligible in 
the tumultuous environment following turnaround. Simi- 
lar hypotheses for A < were previously used in Ref. [2i| . 
It seems reasonable to believe that the truth lies some- 
where between these two hypotheses, perhaps somewhat 
closer to hypothesis B. 



B. Distribution of A using a pocket-based measure 

Before calculating the distribution of A using a scale- 
factor cutoff, we review the standard calculation [24], Hf| 
HI, H3, [H, [29[ . This approach assumes an ensemble of 
equal-size regions with a flat prior distribution of A. The 
regions are allowed to evolve indefinitely, without any 
time cutoff, so in the case of A > the selection factor 
is given by the asymptotic collapse fraction at t — > oo. 
For A < we shall consider anthropic hypotheses A and 
B. This prescription corresponds to using the pocket- 
based measure, in which the ensemble includes spherical 
regions belonging to different pockets and observations 
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FIG. 1: The normalized distribution of A for A > 0, with A 
in units of the observed value, for the pocket-based measure. 
The vertical bar highlights the value we measure, while the 
shaded regions correspond to points more than one and two 
standard deviations from the mean. 



are counted in the entire comoving history of these re- 
gions. The corresponding distribution function is given 
by 



F(M g ,t ^oo) 
P(A) <x { F(Mq, T crunc h — At) 

F{Mq, Tt urn — At) 



for A > 

for A < (A) (23) 
for A < (B), 



where, again, T crunc h = 2ir/3H\ is the proper time of the 
crunch in pockets with A < 0, while T turn = n/SH\. 

We approximate the collapse fraction F using the 
Press-Schechter (PS) formalism [52j], which gives 



F(M g ,t) =erfc 



V2a{M Gl r) 



(24) 



where a(Mc,r) is the root-mean-square fractional den- 
sity contrast SM/M averaged over a comoving scale en- 
closing mass Mq and evaluated at proper time t, while 
5 C is the collapse density threshold. As is further ex- 
plained in appendix [Bl <5 c (t) is determined by considering 
a "top-hat" density perturbation in a flat universe, with 
an arbitrary initial amplitude. S c (t) is then defined as 
the amplitude reached by the linear evolution of an over- 
density of nonrelativistic matter Sp m /p m that has the 
same initial amplitude as a top-hat density perturbation 
that collapses to a singularity in proper time r. S c (t) 
has the constant value of 1.686 in an Einstein-de Sit- 
ter universe (i.e., flat, matter-dominated universe), but 
it evolves with time when A / [54| . We simulate 
this evolution using the fitting functions (|B23[) . which are 
accurate to better than 0.2%. Note, however, that the 
results are not significantly different if one simply uses 
the constant value S c = 1.686. 

Aside from providing the collapse fraction, the PS for- 
malism describes the "mass function," i.e. the distribu- 
tion of halo masses as a function of time. A^-body sim- 
ulations indicate that PS model overestimates the abun- 
dance of halos near the peak of the mass function, while 
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anthropic hypothesis A anthropic hypothesis B 




FIG. 2: The normalized distribution of A, with A in units of the observed value, for the pocket-based measure. The left column 
corresponds to anthropic hypothesis A while the right column corresponds to anthropic hypothesis B. Meanwhile, the top row 
shows P(A) while the bottom row shows P(|A|). The vertical bars highlight the value we measure, while the shaded regions 
correspond to points more than one and two standard deviations from the mean. 



underestimating that of more massive structures [55j |. 
Consequently, other models have been developed (see e.g. 
Refs. [5g|), while others have studied numerical fits to N- 
body results [2i| [57j ■ From each of these approaches, the 
collapse fraction can be obtained by integrating the mass 
function. We have checked that our results are not signif- 
icantly different if we use the fitting formula of Rcf. [29] 
instead of Eq. (J24J). Meanwhile, we prefer Eq. (|24[) to the 
fit of Ref. [29[ because the latter was performed using 
only numerical simulations with A > 0. 

The evolution of the density contrast a is treated lin- 
early, to be consistent with the definition of the collapse 
density threshold S c . Thus we can factorize the behavior 
of <j{Mq,t), writing 

a(M G ,T) = a(M G )G A ( T ), (25) 

where G\(t) is the linear growth function, which is nor- 
malized so that the behavior for small r is given by 
Ga(t) ~ {3H a t/2) 2 / 3 . In appendix E we will give 
exact integral expressions for Ga(t), and also the fit- 
ting formulae (|B1~2"|) and (|B13|) . taken from Ref. [29[, 
that we actually used in our calculations. Note that for 
A > the growth rate Ga(t) always decreases with time 
(Ga(t) < 0), while for A < the growth rate reaches 
a minimum at r w 0.24r crunc h and then starts to accel- 



erate. This accelerating rate of growth is related to the 
increasing rate of matter accretion in collapsed halos af- 
ter turnaround, which we mentioned above in motivating 
the anthropic hypothesis B. 



The prefactor a(Mc) in Eq. (|2"5")l depends on the scale 
Mq at which the density contrast is evaluated. According 
to our anthropic model, Mq should correspond to the 
minimum halo mass for which star formation and heavy 
element retention is efficient. Indeed, the efficiency of 
star formation is seen to show a sharp transition: it falls 
abruptly for halo masses smaller than Mq ~ 2 x IO^Mq, 
where Mq is the solar mass [58[. Peacock [29[ showed 
that the existing data on the evolving stellar density can 
be well described by a Press-Schechter calculation of the 
collapsed density for a single mass scale, with a best fit 
corresponding to ct(Mg, tiooo) ~ 6.74x 10~ 3 , where tiooo 
is the proper time corresponding to a temperature T = 
1000 K. Using cosmological parameters current at the 
time, Peacock found that this perturbation amplitude 
corresponds to an effective galaxy mass of 1.9 x 10 12 Mq. 
Using the more recent WMAP-5 parameters [59j], as is 
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done throughout this paper, 7 we find (using Ref. [6(| and 
the CMBFAST program) that the corresponding effective 
galaxy mass is 1.8 x 10 12 Mq. 

Unless otherwise noted, in this paper we set the pref- 
actor a(Ma) in Eq. J25]) by choosing M G = 1O 12 M . 
Using the WMAP-5 parameters and CMBFAST, wc find 
that at the present cosmic time a(10 12 Mq) w 2.03. This 
corresponds to er(10 12 Mq, nooo) ~ 7.35 x 1CU 3 . 

We are now prepared to display the results, plotting 
P(A) as determined by Eq. ([23| . We first reproduce the 
standard distribution of A, which corresponds to the case 
when A > 0. This is shown in Fig. [TJ We see that the 
value of A that we measure is between one and two stan- 
dard deviations from the mean. Throughout the paper, 
the vertical bars in the plots merely highlight the ob- 
served value of A and do not indicate its experimental 
uncertainty. The quality of the fit depends on the choice 
of scale Mq; in particular, choosing smaller values of Mq 
weakens the fit |28l . |6l| . Note however that the value of 
Mq that we use is already less than that recommended 
by Ref. (H- 

Fig. [2] shows the distribution of A for positive and neg- 
ative values of A. We use At = 5 x 10 9 years, corre- 
sponding roughly to the age of our solar system. The 
left column corresponds to choosing anthropic hypoth- 
esis A while the right column corresponds to anthropic 
hypothesis B. To address the question of whether the 
observed value of |A| lies improbably close to the special 
point A = 0, in the second row we plot the distribu- 
tions for P(|A|). We see that the observed value of A lies 
only a little more than one standard deviation from the 
mean, which is certainly acceptable. (Another measure 
of the "typicality" of our value of A has been studied in 
Ref. Hi]). 



C. Distribution of A using the scale-factor cutoff 

We now turn to the calculation of P(A) using a scale- 
factor cutoff to regulate the diverging volume of the mul- 
tiverse. When we restrict attention to the evolution of 
a small thermalized patch, a cutoff at scale-factor time 
t c corresponds to a proper time cutoff r c , which depends 
on t c and the time at which the patch thermalized, i*. 
Here we take the thermalized patch to be small enough 
that scale-factor time t is essentially constant over hy- 
persurfaces of constant r. Then the various proper and 
scale-factor times are related by 




FIG. 3: The normalized distribution of A for A > 0, with A 
in units of the observed value, for the scale-factor cutoff. The 
vertical bar highlights the value we measure, while the shaded 
regions correspond to points more than one and two standard 
deviations from the mean. 



after which they follow Eqs. (f2"Tj) . Solving for the proper 
time cutoff t c gives 



r c = — H A arcsinh 
3 

for the case A > , and 



— H A ^rcsin 
o 



-c) 



-c) 



(27) 



(28) 



J °H(r) dr = In [o(r )/a(T.)] 



(26) 



for A < 0. The term C is a constant that accounts for 
evolution from time r* to time T m . Note that as t c — t* is 
increased in Eq. (|28[) . r c grows until it reaches the time 
of scale- factor turnaround in the pocket, Tt urn = tt/3Pa, 
after which the expression is ill-defined. Physically, the 
failure of Eq. (|28|) corresponds to when a thermalized 
region reaches turnaround before the scale-factor time 
reaches its cutoff at t c . After turnaround, the scale factor 
decreases; therefore these regions evolve without a cutoff 
all the way up to the time of crunch, r crunc h = 2ir/3H\. 

When counting the number of observers in the various 
pockets using a scale-factor cutoff, one must keep in mind 
the dependence on the thermalized volume V* in Eq. (|20p , 
since in this case U* depends on the cutoff. As stated 
earlier, we assume the rate of thermalization for pockets 
containing universes like ours is independent of A. Thus, 
the total physical volume of all regions that thermalized 
between times £» and f* + dt« is given by Eq. ©, and is 
independent of A. Using Eq. (|20[) to count the number of 
observers in each thermalized patch, and summing over 
all times below the cutoff, we find 



Recall that all of the thermalized regions of interest 
share a common evolution up to the proper time r m , 



P(A) 



F[M G ,T c (t c ,U)- At] e^'dU 



(29) 



7 The relevant values are Q A = 0.742, Q m = 0.258, Q b = 0.044, 
n s = 0.96, h = 0.719, and A^(fc = 0.02 Mpc -1 ) = 2.21 x 10" 9 . 



Note that regions thcrmalizing at a later time i* have a 
greater weight oc e 7 ** . This is an expression of the young- 
ness bias in the scale-factor measure. The A dependence 
of this distribution is implicit in P, which depends on 
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anthropic hypothesis A anthropic hypothesis B 




FIG. 4: The normalized distribution of A, with A in units of the observed value, for the scale-factor cutoff. The left column 
corresponds to anthropic hypothesis A while the right column corresponds to anthropic hypothesis B. Meanwhile, the top row 
shows P(A) while the bottom row shows P(|A|). The vertical bars highlight the value we measure, while the shaded regions 
correspond to points more than one and two standard deviations from the mean. 



<5 c (A,t c - Ar)/cr rms (A, r c - At), and in turn on r c (A), 
which is described below. 

For pockets with A > 0, the cutoff on proper time r c is 
given by Eq. ((27)) . Meanwhile, when A < 0, r c is given by 
Eq. I|28p. when that expression is well-defined. In prac- 
tice, the constant C of Eqs. (f2"7| and (f2"5)) is unimportant, 
since a negligible fraction of structures form before the 
proper time r m . Furthermore, for a reference time r m 
chosen deep in the era of matter domination, the nor- 
malized distribution is independent of r m . As mentioned 
above, for sufficiently large t c — t* Eq. p8|) becomes ill- 
defined, corresponding to the thermalized region reach- 
ing its crunch before the scale-factor cutoff. In this case 
we set r c = r crunc h or t c = T turn , corresponding to the 
anthropic hypothesis A or B described above. 

To compare with previous work, we first display the 
distribution of positive A in Fig. [3] We have set 7 = 3 
and use At = 5 x 10 9 years. Clearly, the scale-factor 
cutoff provides an excellent fit to observation, when at- 
tention is limited to A > 0. Note that the scale- factor- 
cutoff distribution exhibits a much faster fall off at large 
A than the pocket-based distribution in Fig. [TJ The rea- 
son is not difficult to understand. For larger values of A, 
the vacuum energy dominates earlier. The universe then 
begins expanding exponentially, and this quickly triggers 



the scale-factor cutoff. Thus, pockets with larger values 
of A have an earlier cutoff (in terms of the proper time) 
and have less time to evolve observers. This tendency for 
the cutoff to kick in soon after A-domination may help 
to sharpen the anthropic explanation [26|, [62] of the oth- 
erwise mysterious fact that we live so close to this very 
special epoch (matter-A equality) in the history of the 
universe. 

The distribution of A for positive and negative values of 
A is displayed in Fig. 01 using the same parameter values 
as before. We see that the distribution with anthropic 
hypothesis A provides a reasonable fit to observation, 
with the measured value of A appearing just within two 
standard deviations of the mean. Note that the weight 
of this distribution is dominated by negative values of 
A, yet anthropic hypothesis A may not give the most 
accurate accounting of observers in pockets with A < 0. 
Anthropic hypothesis B provides an alternative count of 
the number of observers in regions that crunch before the 
cutoff, and we see that the corresponding distributions 
provide a very good fit to observation. This is the main 
result of this work. 

The above distributions all use At — 5 x 10 9 years and 
Mq = 10 12 Mq. These values are motivated respectively 
by the age of our solar system and by the mass of our 
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FIG. 5: The normalized distribution of A, with A in units of the observed value, for anthropic hypothesis B in the scale-factor 
cutoff. The left panel displays curves for At = 3 (solid), 5 (dashed), and 7 (dotted) xlO 9 years, with Mg = 10 12 Mq, while 
the right panel displays curves for M G = 10 10 M Q (solid), 10 11 Mq (dashed), and 10 12 Mq (dotted), with At = 5 x 10 9 years. 
The vertical bars highlight the value of A that we measure. 



galactic halo, the latter being a good match to an em- 
pirical fit determining the halo mass scale characterizing 
efficient star formation [29(. Yet, to illustrate the depen- 
dence of our main result on Ar and Mq, in Fig. [5] we 
display curves for anthropic hypothesis B, using At = 3, 
5, and 7 x 10 9 years and using M G — 10 10 , 10 11 , and 
10 12 M Q . The distribution varies significantly as a result 
of these changes, but the fit to the observed value of A 
remains good. 



IV. CONCLUSIONS 

To date, several qualitatively distinct measures have 
been proposed to regulate the diverging volume of the 
multiverse. Although theoretical analysis has not pro- 
vided much guidance as to which of these, if any, is 
correct, the various regulating procedures make different 
predictions for the distributions of physical observables. 
Therefore, one can take an empirical approach, compar- 
ing the predictions of various measures to our observa- 
tions, to shed light on what measures are on the right 
track. With this in mind, we have studied some aspects 
of a scale-factor cutoff measure. This measure averages 
over the spacetime volume in a comoving region between 
some initial spacelike hypersurface £ and a final hyper- 
surface of constant time, with time measured in units of 
the local Hubble rate along the comoving geodesies. At 
the end of the calculation, the cutoff on scale-factor time 
is taken to infinity. We shall now summarize what we 
have learned about the scale- factor measure and compare 
its properties to those of other proposed measures. 

The main focus of this paper has been on the proba- 
bility distribution for the cosmological constant A. Al- 
though the statistical distribution of A among states in 
the landscape is assumed to be flat, imposing a scale- 
factor cutoff modulates this distribution to prefer smaller 
values of A. Combined with appropriate anthropic selec- 
tion effects, this gives a distribution of A that is in a 



good fit with observation. We have calculated the dis- 
tribution for positive and negative values of A, as well 
as for the absolute value |A|. For A > 0, we adopted 
the standard assumption that the number of observers is 
proportional to the fraction of matter clustered in halos 
of mass greater than 10 12 M©, and allowed a fixed proper 
time interval Ar = 5 x 10 9 years for the evolution of 
observers in such halos. For A < 0, we considered two 
possible scenarios, which probably bracket the range of 
reasonable possibilities. The first (scenario A) assumes 
that observations can be made all the way to the big 
crunch, so we count all halos formed prior to time Ar 
before the crunch. The second (scenario B) assumes that 
the contracting negative-A phase is hazardous to life, so 
we count only halos that formed at time Ar or earlier 
before the turnaround. 

Our results show that the observed value of A is within 
two standard deviations from the mean for scenario A, 
and within one standard deviation for scenario B. In the 
latter case, the fit is better than that obtained in the 
"standard" calculations [H, [H, US HJ H SI], which 
assume no time cutoff (this is equivalent to choosing a 
pocket-based measure on the multiverse). The causal 
patch measure also selects for smaller values of A provid- 
ing, in the case of positive A, a fit to observation similar 
to that of the scale-factor cutoff [3(J. Note, however, that 
the approach of Ref. [301 ] used an entropy-based anthropic 
weighting (as opposed to the structure-formation-based 
approach used here) and that the distribution of negative 
A has not been studied in this measure. 

We have verified that our results are robust with re- 
spect to changing the parameters Mq and Ar. The 
agreement with the data remains good for Mq varying 
between 10 10 and 10 12 M Q and for Ar varying between 
3 x 10 9 and 7 x 10 9 years. 

We have also shown that the scale-factor cutoff mea- 
sure does not suffer from some of the problems afflicting 
other proposed measures. The most severe of these is the 
"youngness paradox" — the prediction of an extremely 
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youth-dominated distribution of observers — which fol- 
lows from the proper-time cutoff measure. The scale- 
factor cutoff measure, on the other hand, predicts only a 
very mild youngness bias, which is consistent with obser- 
vation. Another problem, which arises in pocket-based 
and "gauge-invariant" measures, is the Q catastrophe, 
where one expects to measure the amplitude of the pri- 
mordial density contrast Q to have an unfavorably large 
or small value. This problem ultimately stems from an 
exponential preference for a large number of e-folds of 
slow-roll inflation in these measures. The scale-factor 
cutoff does not strongly select for more inflation, and thus 
does not suffer from a Q catastrophe. An unattractive 
feature of causal patch and comoving-volumc measures 
is that their predictions are sensitive to the assumptions 
one makes about the initial conditions for the multiverse. 
Meanwhile, the scale-factor cutoff measure is essentially 
independent of the initial state. This property reflects 
the attractor character of eternal inflation: the asymp- 
totic late-time evolution of an eternally inflating universe 
is independent of the starting point. 

As mentioned above, a key features of the scale-factor 
cutoff measure is that, unlike the pocket-based or "gauge- 
invariant" measures, it does not reward large amounts of 
slow-roll inflation. As a result, it allows for the possibility 
of a detectable negative curvature. This issue will be 
discussed in detail in Ref. [49j |. 

With any measure over the multiverse, one must be 
wary that it does not over-predict "Boltzmann brains" 
- observers that pop in and out of existence as a re- 
sult of rare quantum fluctuations 63] . This issue has not 
been addressed here, but our preliminary analysis sug- 
gests that, with some mild assumptions about the land- 
scape, the scale-factor cutoff measure does not have a 
Boltzmann brain problem. We shall return to this issue 
in a separate publication [64j |. 
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APPENDIX A: INDEPENDENCE OF THE 
INITIAL STATE 

In section [TT| we assumed that the landscape is irre- 
ducible, so that any vacuum is accessible through quan- 
tum diffusion or bubble nucleation from any other (de 
Sitter) vacuum. If instead the landscape splits into sev- 



eral disconnected sectors, the scale-factor cutoff can be 

(A) 

used to find the probability distributions in each 

of the sectors (labeled by A). These distributions are 
determined by the dominant eigenstates of the Fokker- 
Planck or master equation, which correspond to the 
largest eigenvalues ja, and are independent of the choice 
of the initial hypersurfaces T,a that are used in imple- 
menting the scale-factor cutoff. But the question still 
remains, how do we compare the probabilities of vacua 
belonging to different sectors? 

Since different sectors are inaccessible from one an- 
other, the probability Pa of being in a given sector must 
depend on the initial state of the universe. For definite- 
ness, we shall assume here that the initial state is de- 
termined by the wave function of the universe, although 
most of the following discussion should apply to any the- 
ory of initial conditions. According to both tunneling 
[6a ] and Hartle-Hawking [66[ proposals for the wave func- 
tion, the universe starts as a 3-sphere S a filled with some 
positive-energy vacuum a. The radius of the 3-sphere is 
r a = H^ 1 , where H a is the de Sitter expansion rate. The 
corresponding nucleation probability is 

PW « exp (± Jy) , (Al) 

where the upper sign is for the Hartle-Hawking and the 
lower is for the tunneling wave function. Once the uni- 
verse has nucleated, it immediately enters de Sitter infla- 
tionary expansion, transitions from a to other vacua, and 
populates the entire sector of the landscape to which the 
vacuum a belongs. We thus have an ensemble of eternally 
inflating universes with initial conditions at 3-surfaces S a 
and the probability distribution P^ cl given by Eq. (IA1[) . 

If the landscape were not disconnected, we could ap- 
ply the scale factor cutoff measure to any single compo- 
nent a of the initial wave function, and the result would 
be the same in all cases. To generalize the scale-factor 
cutoff measure to the disconnected landscape, the most 
straightforward prescription is to apply the scale factor 
cutoff directly to the initial probability ensemble. In that 
case, 

^"AEW'W' (A2) 

Here, 

A/f>(t c )=i4 A )p^e^ (A3) 

is the number of relevant observations in the entire closed 
universe, starting from the hypersurface S a , with a cutoff 
at scale-factor time t c . The are determined by the 
initial volume of the 3-surface S a , and also by the effi- 
ciency with which this initial state couples to the leading 
eigenvector of Eq. ©. In other words, the J\fj(t c ) are 
calculated using S a as the initial hypersurface Sa- Note 
that only the overall normalization of A/'j"' 1 depends on 
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the initial vacuum a; the relative probabilities of different 
vacua in the sector do not. In the limit of t c — > oo, only 
the sectors corresponding to the largest of all dominant 
eigenvalues, 



In 



= max{7 j4 } , 



(A4) 



have a nonzero probability. If there is only one sector 
with this eigenvalue, this selects the sector uniquely. 

Since the issue of initial state dependence is new, one 
might entertain an alternative method of dealing with 
the issue, in which the probability Pa for each sector is 
determined immediately by the initial state, with 



Pa oc £ P t 



(a) 
nucl 



(A5) 



Then one could calculate any probability of interest 
within each sector, using the standard scale factor cutoff 
method, and weight the different sectors by Pa- How- 
ever, although this prescription is well-defined, we would 
advocate the first method that we described as the nat- 
ural extension of the scale factor cutoff measure. First, 
it seems to be more closely related to the description of 
the scale-factor cutoff measure in a connected landscape: 
the only change is to replace the initial state by an en- 
semble of states, determined in principle by one's theory 
of the initial wave function. Second, in a toy theory, 
one could imagine approaching a disconnected landscape 
from a connected one, by gradually decreasing all the 
cross-sector tunneling rates to zero. In that case, the 
limit clearly corresponds to the first description, where 
one sector is selected uniquely if it has the largest domi- 
nant eigenvalue. 

Assuming the first of these prescriptions, the conclu- 
sion is that the probability distribution (|A2|) defined by 
the scale-factor measure is essentially independent of the 
initial distribution (|A1[) . Some dependence on P^ cX sur- 
vives only in a restricted class of models where the land- 
scape splits into a number of sectors with strictly zero 
probability of transitions between them and, in addition, 
where the maximum eigenvalue 7 max is degenerate. Even 
then, this dependence is limited to the relative probabil- 
ity of the sectors characterized by the eigenvalue 7 max - 



APPENDIX B: THE COLLAPSE DENSITY 
THRESHOLD 8 C 

The collapse density threshold 5 C is determined by 
comparing the linearized evolution of matter perturba- 
tions with the nonlinear evolution of a spherical top-hat 
density perturbation, which can be treated as a closed 
FRW universe. The collapse density threshold 5 c (t) is de- 
fined as the amplitude reached by the linear evolution of 
an overdensity 5 = Sp m / p m that has the same initial am- 
plitude as a top-hat density perturbation that collapses 
to a singularity in proper time r. In a matter-dominated 



universe with zero cosmological constant, S c is a con- 
stant; however, it is well known that 5 C depends on the 
collapse time when A is nonzero (see e.g. Refs. (53l |54|). 
In this appendix we first outline the calculation of the 
time evolution of 5 C , then display the results for positive 
and negative A, and finally describe how we apply it in 
our analysis of the collapse fraction F of Eq. (f2"4"|) . 

As suggested by the definition above, both linear and 
nonlinear analyses are involved at different stages of the 
calculation of the collapse density. Arbitrarily small per- 
turbations obey linearized equations of motion, and their 
evolution defines the linear growth function Ga(t): 



5(t) cx G a (t) 



(Bl) 



where Ga(t) is normalized so that the behavior for small 
t is given by G A (r) - (3# a t/2) 2 / 3 , where H A = y^\A\/S. 
The exact nonlinear analysis is used to determine the 
time at which an overdensity with a given initial ampli- 
tude will collapse to a singularity. For simplicity, this is 
worked out for the "top-hat" model, where the overden- 
sity is assumed to be uniform and spherically symmet- 
ric. Such a region is embedded in a flat FRW universe 
containing only non-relativistic matter and cosmological 
constant. 

By Birkhoff's theorem, the evolution of the spherical 
overdensity is equivalent to that of a closed FRW uni- 
verse. The Friedmann equation for a closed FRW uni- 
verse, with scale factor a, may be written as: 



Hi 



sign(A) 



B{k) 



(B2) 



where H = dlna/dr = a/a and B(k) is an arbitrary 
quantity that fixes the normalization of a. We will al- 
ways choose B(0) = 1, so for k = the scale factor is 
normalized in such a way that p m = \p\\ at a = 1 , where 
PA = A/(8tt) is the vacuum energy density. 

Let us first focus our attention on the evolution of a 
linearized density perturbation in a flat FRW universe 
with positive cosmological constant; the case with nega- 
tive cosmological constant proceeds similarly. Consider 
a closed FRW universe obtained by "perturbing" the flat 
universe with a small curvature term 5k. The proper 
time parameter f in such a universe, as a function of the 
scale factor, is given by an expansion with respect to the 
flat background: f (a) = r(a) + 5r(a), where to linear 
order in 5k 



5r(a) 



5k 

2i?A 



dB i 



da' 



(B3) 



(1 + a' 3 ) 3/2 

The scale factor of the closed universe is obtained by 
inverting the function f (a): 

a(r) = a(r) - a{r) 5 T {a{T)) . (B4) 

As mentioned above, the evolution of this closed FRW 
universe also gives the evolution of a small density per- 
turbation. Using p m — (3/87r)_ff^_B(K)/a 3 , one has 



Sp n 



5a dB , , „ dB , . ,„_. 

1— + —(0)=3H5t + —(0), (B5) 
a dK dK 
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where the last equality follows from Eq. (|B4[) . From 
here on, unless noted otherwise, we normalize a so that 
B(k) = 1. It is convenient to introduce the "time" vari- 
able 



Pm 



(B6) 



for both choices of the sign of A. To be consistent with 
Eq. (|B2[) . the solutions for k = are not normalized as 
in Eq. (|2"Tj) . but instead are given by 



a(r) 



sinh 2/3 (|iJ A r) forA>0 



sm 2 / 3 (|iJ A r) 



for A < . 



(B7) 



We can then find the evolution function S(x) from 
Eq. JB5|, using Eq. (JB3J) and also Eq. jB2| with k = 0: 



dy 



3( x ) = 7^ K \/l + _ ; , , , •■■■■> 

2 V X J y 1 /6(l + y)3/2 

= -SkG+(x), 
5 

where the linear growth function (for A > 0), 

dy 



.„)3/a 



(B8) 



(B9) 



is normalized so that the behavior for small x is given by 
G+{x) ~ x 1 / 3 = a ~ (3# A r/2) 2 / 3 . 

In the A < case, the calculation proceeds along the 
same steps as before and the formula (|B8|) is indeed valid 
also for negative A, after replacing the growth function 
with G~(x). This function now has two branches Gi(x) 



and G^x), corresponding to the expanding and con- 
tracting phases of the universe, respectively. The first 
branch of the growth function introduces no new compli- 
cations, and is found to be 



cr(x) 



dy 



yl/e(l-y)3/2 



(BIO) 



For the second branch, the integration is first performed 
over the whole history of the universe, from x = to 
x = 1 and back to x = 0, and then one integrates back to 
the value of interest x. There is a complication, however, 
because for this case the denominator in Eq. (|B3[) is (1 — 
a' 3 ) 3 / 2 , so the integral diverges when the upper limit is 
equal to 1. The cause of the problem is that for 6 k ^ 
0, a max is no longer equal to 1. A simple cure is to 
choose B(k) = 1 + k for this case, which ensures that 
flmax = 1 for any k, and which correspondingly provides 
an additional term in Eq. (|B3p which causes the integral 
to converge. After some manipulation of the integrals, 
the result can be written as 



dy 



yl/6 (1 _ y) 3/2 



The time dependence of the linear growth functions can 
be made explicit by expressing i as a function of r, 
through Eqs. |B6| and l[B"7|) . 

In practice, we carry out our calculations using fitting 
functions for the growth functions, which were devised 
by Peacock [2^|, and which are accurate to better than 
0.1%. These give 



Gl(r) 
G7(t) 



tanh 2/3 (|i7 A r) [l - tanh 127 '(|-ff A r) 
(|^ A r) 2/3 [l + 0.37(r/r crunch ) 2 - 18 



0.82 



1.4375 



-2/3 



l-cosh- 4/3 (§H A r) 



1 (^"/^"crunch) 



(B12) 
(B13) 



r 



for the cases A > and A < 0, respectively, where the 
latter fitting formula is valid for both branches. 

We are now prepared to set the calculation of 5 C . Since 
the universe in Eq. (|B2j) can be viewed as a "perturba- 
tion" over a flat universe with 5k — k, the time evolution 
of the overdensity is described in general by 



5{t) = -kG a (t). 
5 



(B14) 



The quantity (3/5) k a quantifies the size of the initial 
inhomogeneity. 

In order to find the time at which the spherical over- 
density collapses, it is convenient to determine the time 
of turnaround Tt urn , corresponding to when H = 0. The 



time of collapse is then given by 2 Tturn- The turnaround 
time is obtained by integrating Eq. (|B2[) , choosing B = 1 : 



H A T tmn (K) = 



aturn(n) 



fa da 



o y sign(A) a 3 — k a + 1 



(B15) 



where the scale factor at turnaround atum corresponds 
to the smallest positive solution of 



sign(A) a 3 urn - k a turn + 1 = 0. 



(B16) 



For positive A, the universe will collapse only if k > 
K min = 3/2 2 / 3 ; for negative A, perturbations that col- 
lapse before the universe has crunched have k > 0. 
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2.0 - 




FIG. 6: The collapse density thresholds <5+ (for A > 0) and 
<5~ (for A < 0), as functions of time. The solid curves rep- 
resent numerical evaluations of <5*, while the dashed curves 
correspond to the fitting functions in Eq. (|B23[) . Note that 
5f decreases with time, while 8~ increases with time. 



The numerical evaluation of the integral in Eq. (|B15|) 
allows one the extract the function Tt urn (K), which can 
be inverted to give K turn (r), expressing the value of k 
that leads to turnaround at time r. Finally, the collapse 
density threshold as a function of the time of collapse is 
read from Eq. (|B14|) : 



5 c {t) = -K t um(r/2) G a (t) 
5 



(B17) 



In the limits of small and large collapse times the above 
procedure can be carried out analytically to find the lim- 
iting values of 5 C . Let us consider first the large-time 
regime, corresponding to small k. In the case A > 0, the 
smallest k allowed is K mi „: therefore 



5+M 



.G+M- 1.629. 



(B18) 



where G^(oo) = G+(oo) = 5r(2/3)r(5/6)/(3V^) - 
1.437. The case A < is a little more complicated. The 
collapse time cannot exceed r crunc h = 2tt/3H\, corre- 
sponding to K = 0. At small k, the integral in Eq. (|B15|) 
is expanded to give 



H\T tu 



7 -^A ^crunch 



s r (§) 



(B19) 



On the other hand, the growth function G (r) in the 
neighborhood of r crunc h behaves as 



( T ~ ^crunch) 



io r(|) 



i 



3 (i) (1 - T/T crU nch) 



(B20) 



After using Eqs. (|B19|1 and (|B20[1 in the general formula 
(IB17|) . we simply get 



&c (^crunch) 2 . 



(B21) 



In the opposite regime H\t <C 1, corresponding to 
large k, the growth functions are G^(t) ~ a(r) ~ 
(3H a t/2) 2 / 3 . The integral (|B15| can be analytically 
solved in this limit: H\T tUTn (n) = 7t/(2k 3 / 2 ). Combining 
these results leads to 



2/3 



1.686, 



(B22) 



which is also the constant value of 6 C in a A = universe. 

The time dependence of S c is displayed in Fig. [SJ for 
both positive and negative values of A. We also display 
the following simple fitting functions, 



-0.28HiT J 



S+( T ) = 1.629 + 0.057 e" 
5~{t) = 1.686 + 0.165 ( T ) 

\ ^crunch / 



2.5 



0.149 



li 



\ ^crunch / 

(B23) 



which are accurate to better than 0.2%. Although we 
choose to include the effect of the time evolution of S c , 
our results are not significantly changed by treating S c 
as a constant. This is easy to understand. First of all, 
5£ varies by only about 3%. The evolution of S ~ is more 
significant, about 15%, and most of this happens at very 
late times. But our anthropic weight in Eq. (|20p never 
samples 8 ~ within a time At of r crunc h. 

Finally, we point out that the appearance of Ga(t) 
in this discussion is not needed for the calculation, and 
appears here primarily to make contact with other work. 
From Eq. (pM)) one sees that the collapse fraction depends 
only on the ratio of 6 c (t)/<j(Mg, t), which from Eqs. (T2"5")) 
and (|B17P can be seen to equal (3/5)Kt um (r/2)/0'(MG)- 
Expressed in this way, Eq. (|24|) becomes fairly transpar- 
ent. Since k is a measure of the amplitude of an initial 
perturbation, Eq. (pM]) is saying that the collapse fraction 
at time r depends precisely on the magnitude required 
for an initial top-hat perturbation to collapse by time 
t. In more detail, Eq. flH]) is predicated on a Gaussian 
distribution of initial fluctuations, where the complemen- 
tary error function erfc(a;) is the integral of a Gaussian. 
The collapsed fraction at time r is given by the prob- 
ability, in this Gaussian approximation, for the initial 
fluctuations to exceed the magnitude needed for collapse 
at time r. From a practical point of view, the use of 
Ga(t) in the discussion of the collapse fraction can be a 
helpful simplification if one uses the approximation that 
5 C « const. We have not used this approximation, but as 
described above, our results would not be much different 
if we had. We have maintained the discussion in terms of 
Ga(t) to clarify the relationship between our work and 
this approximation. 
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